Many-Body Effects on Nonadiabatic Feshbach Conversion in Bosonic Systems 
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We investigate the dynamics of converting cold bosonic atoms to molecules when an external 
magnetic field is swept across a Feshbach resonance. Our analysis relies on a quantum microscopic 
model that accounts for many-body effects in the association process. We show that the picture 
of two-body molecular production depicted by the Landau-Zener model is significantly altered due 
to many-body effects. In the nonadiabatic regime, we derive an analytic expression for molecular 
conversion efficiency that explains the discrepancy between the prediction of the Landau-Zener 
formula and the experimental data[Hodby et al., Phys. Rev. Lett. 94, 120402 (2005)]. Our theory 
is further extended to the formation of heteronuclear diatomic molecules and gives some interesting 
predictions. 
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I. INTRODUCTION 

The production of ultracold diatomic molecules is 
an exciting area of research with important applica- 
tions ranging from the search for the permanent elec- 
tric dipole moment [1] to BCS-BEC (Bose-Einstein con- 
densate) crossover physics 0- A widely used produc- 
tion technique involves the association of ultracold atoms 
into very weakly bound diatomic molecules by applying 
a time varying magnetic field in the vicinity of a Fes- 
hbach resonance [Oil- The underlying conversion dy- 
namics are usually described by the Landau-Zener (LZ) 
model[5fl. In this model, the Feshbach molecular produc- 
tion is discussed under a two-body configuration where 
a single pair of atoms is converted to a molecule at an 
avoided-crossing between atomic energy level and molec- 
ular energy level while the molecular energy is lifted by 
an applied linearly sweeping magnetic field. Thus, the 
molecular production efficiency is expected to be an ex- 
ponential Landau-Zener type(y,0|- However, recent ex- 
perimental data on 85 Rb by the JILA group showed a 
large discrepancy from the Landau-Zener formula: The 
value of the LZ parameter extracted from the data is 8 
times larger than the prediction of LZ theory Q. The 
experiment was performed under unusually low densities 
of the atom cloud ( ~ 10 11 /cm 3 ) and the data was mea- 
sured in the nonadiabatic regime so that the inverse ramp 
rate was less than 100/is/G. Therefore, two- and three- 
body atomic decay and collisional molecular decay rates 
are negligible and do not affect the measurement. This 
puzzle remains unresolved and challenges our knowledge 
of the big issue of Feshbach molecular formation. 

In this paper, using a many-body two-channel micro- 
scopic Hamiltonian, we investigate the dynamics of Fes- 
hbach molecular formation in bosonic systems such as 
85 Rb. We show that many-body effects alter the LZ 
picture of two-body molecular production through dra- 
matically distorting the energy levels near the Feshbach 



resonance. With the help of a mean-field classical Hamil- 
tonian, we derive an analytic expression for the conver- 
sion efficiency in the nonadiabatic regime. Our theory 
agrees with experimental data. Our theory thus is ex- 
tended to the Feshbach formation of heteronuclear di- 
atomic molecules such as 85 Rb- 87 Rb and predicts that 
many-body effects are more significant there. 

Our paper is organized as follows. In Sec. II, we present 
our model. In Sec. Ill, many-body effects on the conver- 
sion dynamics are addressed and an analytic expression 
for conversion efficiency is derived. In Sec. IV, we apply 
our theory to explain the experimental data. In Sec.V, we 
extend our theory to the heteronuclear molecules. The 
final section is our conclusion. 



II. MODEL 

Considering the experimental condition that the den- 
sities of the atom cloud is unusually low and the two- and 
three-body atomic decay and collisional molecular decay 
rates are negligible, we exploit the following two-channel 
model to describe the dynamics of converting atoms to 
molecules in the bosonic system, 



H = e a a}a + e b (t)Pb + -jL [a'a^b + Paa 
\/V \ 



(1) 



where a (a + ) and b(b + ) are Bose annihilation (creation) 
operators of atoms and molecules, respectively, g = 
\J 4:Trh 2 ab g ABfi co /m is the atom-molecule coupling due 
to the Feshbach resonance, m is the mass of a bosonic 
atom, dbg is background scattering length, AB is the 
width of the resonance, and fi co is the difference in the 
magnetic moment between the closed channel and open 
channel state. Here, we introduce parameter V to denote 
the volume of trapped particles and therefore n = N/V is 
the mean density of initial bosonic atoms. The external 
magnetic field is linearly swept B(t) = Bt and crosses the 
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FIG. 1: Schematic plot of energy levels of atomic mode and 
molecular mode. Our single-mode approximation is valid only 
when the energy distribution of the thermal particles is much 
smaller than the effective width of the Feshbach resonance, 
i.e., fcsT << ^/g 2 n. Moreover, it is required that the time 
for the external magnetic field to sweep across the Feshbach 
resonance is smaller than the "dephasing" time Td = jrpf- 
For a detailed discussion refer to the text. 



Feshbach resonance at Bq. The molecular energy under 
the field is £f,(t) = fj, co (B(t) — B ). The total number of 
particles N = a T a + 2tfb is a conserved constant. 

The above single-mode model is an approximate de- 
scription of the clouds of noncondensed atoms and 
molecules, and is only valid when the energy distribu- 
tion of the thermal particles (characterized by fc^T, ks is 
the Boltzman constant and T is the temperature) is much 
smaller than the effective Feshbach resonance width g^/n. 
In such cases, each 'energy band' of the thermal parti- 
cles can be approximately denoted by one energy level, as 
schematically plotted by FigfTJ Initially, the particles on 
one level have a definite phase and the phase difference 
between two levels is well defined. However, as the mag- 
netic field sweeps across the Feshbach resonance from B ln 
to B en d at a rate of B, particles will acquire additional 
phases that are proportional to their individual energy 
and sweeping time (B en d — Bi n )/B. The varied particles 
in one level could acquire different phases because they 
have different energies. The validity of our single-mode 
approximation requires that the above mismatch in phase 
or 'dephasing' is at least smaller than 2ir, which defines 
a 'dephasing time' = When the time taken by 

the external magnetic field to sweep across the Feshbach 
resonance is smaller than the above 'dephasing' time, i.e., 
{Bend — Bi n )/B < j^ij, , the above dephasing effects can 
be ignored. The above analysis sets up a lower bound 
on the sweeping rate. So, in the following discussion, we 
focus on the fast-swept or nonadiabatic regime in which 
the above condition is satisfied. 

Using the Fock states as a basis, the Schrodinger equa- 



FIG. 2: Energy levels versus the scaled external magnetic 
fields, a) Homonuclear two-body case; b) Homonuclear many- 
body case (TV = oo), the dashed line represents an ad- 
ditional unstable eigenstate; c-d) Heteronuclear many-body 
cases (TV = oo), where parameter d denotes initial population 
imbalance between two species. 

tion is written as 

i±W) = H\1>)> (2) 

where |V> = Eflo ^\2j, N/2 - j), \2j,N/2 - j) = 
- 1 (tftf) 3 Itf) ' |0) (j = 0,...,/V/2) 

are Fock states, and Cj is the probability amplitudes on 
the corresponding Fock state, respectively. The normal- 
ization condition is Ej \ c j\ 2 = 1- 

For the simplest case of /V = 2, the above Schrodinger 
equation reduces to the following two-level system of 
Landau-Zener type, 

*l (£)-(?/» £)(£)■ (3) 

Where |crj| 2 and |ci| 2 denote the population of molecules 
and atoms, respectively. For the two-level system, the 
energy bias between two levels is 7 = (2e a — e&) and the 
coupling strength is given by v — 2g^/n. Initially, all 
particles populate in the lower level of the atomic state, 
i.e., c = 0,Ci = 1. When the external magnetic field is 
linearly swept across the Feshbach resonance at 7 ~ 0, 
a fraction of atoms will be converted to molecules at the 
avoided-crossing of energy levels. The conversion effi- 
ciency as a function of the sweeping rate (i.e., 7 = /J, C oB) 
and coupling strength, takes the formQ, 

T lz = l- exp -— = 1 - exp L .4 

2?i7 V m\B\ J 

The above is the two-body molecular production picture 
and is consistent with the result from the coupled-channel 
scattering calculation in Ref.[7|. 
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FIG. 3: Feshbach molecular production efficiency versus the 
scaled inverse sweep rates. As we increase the number N of 
the particles, the many-body result converges to the mean- 
field curve. 



Mathematically, ignoring a total phase, the dynamics 
of Eq. ([3]) ar e eq uivalent to the following simple classical 
Hamiltonian LLCl FlOI , 



Hi. 



-"f/hs + v jfi\pY 



! cos 9. 



(5) 



Where the canonical conjugate variables are the popu- 
lation difference s = | Co | 2 — |ci| and the relative phase 
9 = arg Co — arg ci . The dynamics are governed by the 
canonical equations of 8 = , s = — gj- . The fixed 

points satisfying s = 0,9 = correspond to the ex- 
ternum of system energy. These classical fixed points 
correspond to the eigenstates of quantum equations (|3|) 
and their energies (corresponding to the eigenvalues of 
quantum eigenstates) are calculated and plotted against 
the energy bias parameter 7 in Fig(2^i. It exhibits a typ- 
ical LZ avoided-crossing configuration. Initially, all par- 
ticles populate in the atomic state of sq = — 1 at the left 
end of the lower level. When the external field passes 
through the Feshbach resonance of width v/K at 7 = 0, a 
fraction of atoms are converted to molecules at the right 
end of the lower level, leading to a variation in the pop- 
ulation variable, i.e., 



s f = 2T lz - 1 = l-2exp(-— ). (6) 

z«7 

As we go beyond the above two-body treatment to con- 
sider the many-body effects, the structure of the energy 
levels will change dramatically and the above LZ formula 
of the conversion efficiency will be altered due to many- 
body effects. 



III. MANY-BODY EFFECTS IN FORMING 
HOMONUCLEAR FESHBACH MOLECULES 



To include many-body effects, we need to solve full 
+ 1 dimensional quantum equations Using the 



basis of Fock states, the Schrodinger equation is rewritten 
as 

i^ = J2 H ^k,(j,k = 0,l,-,N/2) (7) 

k 

where the Hamiltonian 
2j,N/2 - j\H\2k,N/2 - 
3T, for 3 ^ k , Hjk = 



matrix elements 
- k >. For j ■■ 



e H jk 
k, H 



=< 



except H 



H 



33 



^/(j + l)(2j + l)(N/2-j)/2Nv. 

The above differential equations do not have explicit 
analytic solutions. We thus solve them numerically using 
the 4-5th order Runge-Kutta algorithm with an adaptive 
time-step. Our result is presented in FigJ31 which shows 
that the molecular conversion is altered due to many- 
body effects. Interestingly, there exists a crossing point 
between the two-body conversion curve and the many- 
body conversion curve. As the scaled inverse sweep rate is 
below this crossing point, the many-body effects enhance 
the molecular conversion efficiency, while as the scaled 
inverse sweep rate is above the crossing point, the many- 
body effects suppress the molecular conversion efficiency. 
The location of the crossing point is dependent on the 
total particle number N and shifts to the right as N 
increases. For N = 6, the crossing point corresponds to 
v 2 /tvy = 0.7. It shifts to one as N — 00. 

Below, with the help of angular momentum operators, 
we deduce an analytic expression for the atom-molecule 
conversion efficiency under the mean field approximation. 

The ang ular momentum operators are introduced as 
follows Ha, 



a* cvb + b'aa 



TV 3 / 2 



aJcJb — b'aa 



L z = 



v/2 

2tib-tfa 
N ' 



(8) 
(9) 
(10) 



The operator L z denotes the atom-molecule popula- 
tion imbalance, and L x ,L y describe the coherence of 
atoms and molecules. They compose a generalized Bloch 
representation The commutators between the opera- 



tors are 



L 2 , hr. 



a f 

N V 



L Z , Ly 



4l T 



L X 1 Ly 



jf (l - L^j (l + 3L Z J + L x ,L yi L z do not span 



SU(2) because the commutator 



L X ,L 



yields a 



quadratic polynomial in L z . The generalized Bloch 

surface is determined by the conserved relationship 

2 



{L x y + (L v Y = \[l + L z ) [l-l 

-£tL z . The Hamiltonian (fT]) becomes H = —^.jL z + 

^^-vL x ^r^. The Heisenberg equations are ih-^Lj = 

Lj,H ,j = x,y,z. 

In the mean field limit where N — > 00, all the above 
commutators vanish. Therefore, it is appropriate to 
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replace L x , L y and L z by their expected values u,w, 
and s, respectively. Noting the constraint u 2 + w 2 — 
|(s — l) 2 (s + 1) and introducing the conjugate angular 
variable 8 = arctan(w/it) denoting the relative phase be- 
tween atoms and molecules, the Heisenberg equations can 
be replaced by a classical Hamiltonian of the form 



H r , 



-7/fts + v/hy/(l - s 2 )(l - s) cos9, (11) 

To understand the dynamics, we first look at the fixed 
points s = 9 = 0. The energies for these fixed points 
make up energy levels of the system, as shown in Fig|2jD. 
The structure of these energy levels changes dramatically 
compared to the two-body case. We observe: (i) There 
are two fixed points when \~f/v\ is large enough: one for 
the bosonic molecule (BM) and the other for the bosonic 
atom (BA). (it) When \j/v\ < y2j there is an addi- 
tional fixed point with s = 1. However, this fixed point 
is a saddle point corresponding to dynamically unstable 
quantum states[15]. Eq. (|1 1|) and the energy spectra in 
FiglSb are the same as those obtained for the two-mode 
atom-molecule Fermi system[l3|, [l6| except that the sign 
of the magnetic field is reversed. 

Compared to Hamiltonian ([5]), the coupling term in 
many-body Hamiltonian (|11|) is renormalized by a fac- 
tor vl — s. So, the Fehsbach resonance width that is 
proportional to the coupling either broadens or shrinks 
depending on the factor. For the fast sweep case, s should 
be not far from its initial value —1, therefore, the reso- 
nance width broadens and we expect that many-body ef- 
fects enhance the atom- molecule conversion. In contrast, 
for the slow sweep case, s may take a value close to 1, 
therefore, the resonance width shrinks. We then expect 
that the many-body effects suppress the atom-molecule 
conversion compared to the two-body Landau- Zener for- 
mula. The above analysis reveals the mechanism behind 
the crossing phenomenon exhibited in Fig[3J 

To derive an approximate analytic expression for the 
conversion efficiency, we introduce an effective coupling 
v eff as, 



"eff 



(12) 



Where s* can be approximately taken as the average 
between initial value sq = — 1 and the final value s/, 
i.e., s* = (—1 + Sf)/2. Using the relation T m = 2 < 
6+6 > /N = (1 + s f )/2 and formula (4), we obtain a 
self-consistent formula for the many-body conversion ef- 
ficiency r m , 

irv 2 (2-T m ), 



r, 



i 



exp( 



2hj 



-)■ 



(13) 



The above self-consistent equation for the conversion ef- 
ficiency T m can be readily solved using the iteration 
method. The result is presented in FigJ3] We also nu- 
merically solve the mean field equations using the Runge- 
Kutta step-adaptive algorithm. They are in good agree- 
ment, especially in the nonadiabatic regime of fast sweep 
rates (see Fig(3J) . For the slow sweep case, the above 
formula overestimates molecular conversion slightly. 
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FIG. 4: Molecular conversion efficiency versus the inverse 
sweeping rates. 



IV. COMPARISON WITH EXPERIMENT OF 



RB 



Now wc apply our theory to 85 Rb experiment by the 
JILA group Q . The atoms are held in a purely magnetic 
"baseball" trap. For efficient evaporation, the bias field 
is held at 162 G, where the scattering length is positive. 
For slow magnetic field ramps, Rb2 molecules are pro- 
duced only when the field is ramped upward through the 
resonance, which is located at 155G. Hence, the first step 
in molecule production is to rapidly jump the magnetic 
field from 162G to 147. 5G. They then sweep the field back 
up to 162G at a chosen linear rate, producing molecules 
as they pass through the Feshbach resonance. The ini- 
tial conditions of the atomic cloud are N — 87000 and 
n = 1.3 x 10 11 cm -3 . The Feshbach resonance parameters 
are db g = — 443an,A£> = 10.71G, and \i co = — 2.33^s, 
where an and [Ib are the Bohr radius and Bohr magne- 
ton, respectively. The thermal cloud of the particles is at 
temperature T = 40nK. 

In the first part of their experiment, they measured 
the molecular conversion efficiency as a function of the 
inverse ramp rate. A typical data set is shown in Fig|?] 
The data are mainly divided into three regimes, i.e., the 
linear increase regime where the inverse ramp rates are 
less than 100 /is/G, the saturation regime where the 
inverse ramp rates are larger than 200 /is/G, and the 
transition regime in between. The single-mode approx- 
imation exploited in our theory requires that the reso- 
nance width is much larger than the energy distribution 
of the particles. In this experiment, we see that the ratio 
g^/n/k B T is around 20 at T = 40 nK. While in the sat- 
uration regime, the molecular conversion rates are found 
to saturate around 37%. From our "dephasing" criterion 
discussed in Sec. II, i.e., (B en d — B in )/B < jj^Sp and that 

B m = 147.5G, B e 
So, 



= 162G, we have l/B < 82/is/G. 
in our discussion, we only focus on the first regime. 
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FIG. 5: The ratios of mean density over 1/e ramp rate, with 
respect to mean density. Our theory shows a good agree- 
ment with the experiment for four low density points but is 
obviously larger than the final point. At the high density, 
the cloud experienced significant heating during the ramps 
across the resonance, hence the density of the final point has 
significant uncertainty^. 



In Fig 01 we plot the results from our many-body the- 
ory, which show good agreement with the experimental 
data in the linear regime. As a comparison, we also plot 
the result from the Landau-Zener formula, which shows 
a pronounced deviation from the experimental data. 

In the second part of the experiment, to compare 
with Landau-Zener theory quantitatively, the JILA group 
measured the ratio between mean density and 1 /e ramp 
rate as a function of mean density. They found that 
the Landau-Zener parameter predicted from the two- 
body theory is roughly 1 /8 of the value extracted from 
the experimental data. They use the formula N mo i = 
N max {l - e -omABa bg /B^ to nt the experimental data 
on molecular conversion, where N max is the asymptotic 
number of molecules created for a very slow ramp,B is 
the magnetic field sweeping rate, and a is a fitting pa- 
rameter. Slz — anABabg/B is the Landau-Zener pa- 
rameter. The saturation data in Fig|4] indicate that 
N max /N = 37%. 
as that at B 1 / e , S 



The 1/e ramp rate B 1 / e is defined 
lz = 1 and N mo i/N max = 63%. 
It was then claimed that the data support a constant 
value for n/Bi/ e (see FigJS]). The value for a, extracted 
from the experimental data, is 4.5 x 10~ 7 m 2 s~ 1 . How- 
ever, the two-body Landau-Zener formula (4) predicts 
a = 8ir 2 h/m = 5.9 x 10 -8 m 2 s -1 , roughly 1/8 of the 
experimental data. 

Now we apply our many-body theory to resolve this 
puzzle. At -Bi/ e , the molecular conversion efficiency 
is N mol /N = 37% x 63% = 23%. In the nonadia- 
batic regime, our many-body formula Q13p is simplified as 

r m ~ gg?gte£gl . Substituting r m = 23%, B = B 1/e 
into the above formula, we have n/Bi/ e — w ^ ) h\aT&B\ = 



105 x 10 11 cm 3 fisG 1 , which is good agreement with the 
experimental data of the fourth scatter in FigJS] 

To compare with two-body LZ formula Eq.(4), we see 
that, the many-body effects change the 1/e rate in the 
non-adiabatic regime by a factor of 2. The above analysis 
uncovers the physics behind the 1/8 deviation. The fac- 
tor 1/8 is the product of following three factors: 0.37 is 
from the maximum conversion rate, 0.63 is from the def- 
inition of the 1/e ramp rate, and 1/2 comes from many- 
body effects. . 

Our calculations are extended to the cases of varied 
spatial densities. As mentioned above, in Ref. 8], the for- 
mula N mo l = N max (l ~ e -"»ABa b3 /B) j g uged tQ fit ex _ 

pcrimental data on molecular conversion. Accordingly, 
the 1/e ramp rate B\u corresponds to N mo i/N max = 
1 — 1/e = 63%. Our many-body theory Eq.(13) predicts 
that n/B 1/e = le7T ° h ^" g \ B \ Because the maxi- 

mum molecular conversion efficiency (i.e., N max /N) is 
a function ofpeak phase space density as revealed in 
Fig. 2 in Ref. 3 and the spatial density is proportional 
to peak phase space density at the fixed temperature, we 
claim that n/Bi/ e is spatial density n dependent through 
N m ax/N . The N max /N as a function of density is read 
out from Fig. 2 in Ref.[8[. Thus, our theoretical curve is 
plotted against the experimental data in Figj5] It shows 
good agreement with the experiment for four low density 
points but is obviously larger than the final point. At 
high density, the cloud experienced significant heating 
during the ramps across the resonance, hence the den- 
sity of the final point has significant uncertainty (i.e., see 
the caption of Fig.l of Q). The result from two-body 
LZ theory is also presented in FigOfor comparison. It is 
twice as large as that of many-body theory, and obviously 
deviates from the experimental data. 



V. MANY-BODY EFFECTS IN FORMING 
HETERONUCLEAR FESHBACH MOLECULES 



In the above discussion, we investigated the dynamics 
of Feshbach converting single atomic species to homonul- 
cear diatomic molecule. Actually, the Feshbach reso- 
nance technique has been used to produce heteronuclear 
molecules from two or more species of atoms p~7| . These 
ultracold heteronuclear molecules in low-lying vibrational 
states are of particular interest since they could be a per- 
manent dipole moment due to the unequal distribution of 
electrons. Other proposals for using the polar molecules 
include quantum computation[18j and testing fundamen- 
tal symmetry [l9j]. 

In this section, we extend our discussion to the two- 
species atom case and show that the heteronuclear molec- 
ular production efficiency is more significantly altered 
due to many-body effects. The many-body three-channel 
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Hamiltonian for the heteronuclear system reads [2C 



H 



e a 2a 2 d 2 



-y= (ji\a\b + 6^01(22) . 



(14) 



where a\ (a\), a 2 (aj) are annihilation (creation) op- 
erators of the heteronuclears atoms and b (b') are 
annihilation (creation) operators of molecules, g — 
y^2irh 2 abgAB^, co /'m / is the atom-molecule coupling 
strength, and ml = m\m 2 l (m\+rn 2 ) is the reduced mass 
of two atom scattering. The total number of particles 
N = N a i + N a2 + Nb = a\ai + d\d 2 + 2bflb is a conserved 
constant. The density n = N/V. 

Using the Fock states as a basis, the Schrodinger equa- 
tion is written as 



(15) 



where 



y/j\j\(N/2-jy. 



E"l 2 oC 3 \3,3,N/2 - j), \j,j,N/2 j) = 
■ ' llA ■ |0) (i = 0,..., N/2) are 



Fock states, and Cj is the probability amplitudes on the 
corresponding Fock state, respectively. The normaliza- 
tion condition is that £\ \cj\ 2 = 1. 



For N — 2, the Schrodinger equation reduces to 



ih 



dt 



e b _ v/2V2 
v/2y/2 e al + e a2 



(16) 



Where |co| 2 and |ci| 2 denote the population of molecules 
and atoms, respectively. The energy bias 7 = e a i + e Q 2 — 
t b . Then, two-body molecular production efficiency is 1 — 
exp(— jjp-). Comparing the above expression with Eq.(j4]), 
a 1/2 factor emerges in the exponent. This is due to 
the distinguishability between two atomic species that 
decreases the effective density of each atomic species. 

For the two-species case, the number of particles in 
each species may not be identical. Therefore, we intro- 
duce a parameter d to denote the population imbalance 
between the two species, i.e., d = (iV al — N a2 )/N as- 
suming that N a i > N a2 . Our concern is the conversion 
efficiency of type-2 atoms, i.e., Td — 2N b /{l — d)N when 
the magnetic field is swept across the resonance. We will 
show that d is an important parameter in forming the 
heteronuclear molecule. For the larger population im- 
balance, the heteronuclear molecule production is more 
significantly altered due to the many body effect. 

For the heteronuclear system, the Bloch space 
is expanded by following three operators L x = 



2V2- 



AT3/2 



and L, — 



' jy3/2 7 Zy Zl = 

2 ^ k ° iai — < h22. The commutators between the operators 



L Z i L x 

L*X1 Ly 



I 

N 



(!-!,)( 



1 + 3L ; 



Ad 2 



Ai 



l + L t 



(17) 
(18) 

(19) 



The generalized Bloch surface is determined by the con- 
served relationship 



{L X Y + (L V Y = [1 + L Z + 



N 



1-L t 



Ad 1 



(20) 



Hamiltonian (11411 becomes H 



£ 7 Z Z + "vL 



The 



Lj,H 



4V2 
,3 = x,y,z. 
it is appro- 



- — E 
4 

Heisenberg equations are ih-^Lj = 

In the mean field limit where N 
priate to replace L x , L v , and L z by their expected 
values u,w and s, respectively. Noting the constraint 
u 2 + w 2 = (1 + s) ((1 — s) 2 — Ad 2 ) and introducing the 
conjugate angular variable 6 = arctan(u>/u) denoting the 
relative phase between atoms and molecules, the Heisen- 
berg equations can be replaced by a classical Hamiltonian 
of the form 

K = ~\s + JjU y/(l + s) [(1 - s) 2 - Ad 2 } cos9, (21) 



and the canonical equations, 



dd/dt 
ds/dt 



cos(fi,2) 



7 v (l-s)(l + 3s)+4d 2 

~h ~ 2V2h ^/JTTs) [(1 - s) 2 - aW] 

JL_y/(l + s) [(l- S )2-4^]sin(0), (23) 



The fixed points of the above system have been ob- 
tained by setting s = 6 = 0. The energies for these 
fixed points make up energy levels of the system, as 
shown in Fig[2j:,d. The structure of these energy lev- 
els changes dramatically compared to the homonuclear 
case. There are always two fixed points corresponding 
to two branches of energy levels. Moreover, for a large 
population imbalance between two species, for example 
d = 0.9, the energy levels tend to parallel each other 
(see. Fig|3JI). Thus, the level space remains almost con- 
stant and is slightly dependent on the external field. We 
therefore expect that molecular efficiency is very large 
in this case. Compared to Hamiltonian ([5]), the cou- 
pling term in many-body Hamiltonian (|21[) is renormal- 

ized by a factor of 



^ 2(l-s) d • ^° derive an approx- 
imate expression, we use the effective coupling v e ff as, 



V eff 



(l-s-) 2 -4ri 2 
2(W) 



, where s* can be approximately 



N 



taken as the average between the initial value sq = — 1 
and the final value Sf, i.e., s* = (—1 + s/)/2. Using the 
relation T d = 2N b /N(l - d) = (1 + s/)/2(l - d) and 
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FIG. 6: Feshbach heteronuclear molecule production effi- 
ciency versus the scaled inverse sweep rate. The solid curves 
are our analytical result. For details refer to the text. 



formula ([6]), we obtain a self-consistent formula for the 
many-body conversion efficiency T^, 

(1 - d)T d ~ 

1 - exp ?_ W(2( 1 -d»)-r,(l- < i)(l + «P)) \ 

V Ah i / 

The above self-consistent equation for the conversion 
efficiency Yd can be readily solved using the iteration 
method. The result is presented in Fig|BJ We also numer- 
ically solve the mean field equations for comparison using 
the Runge-Kutta step-adaptive algorithm. The agree- 
ment is good in the nonadiabatic regime where the con- 
version rate is less than 0.6 (see FigJB]). For d = 0.9, the 
regime corresponds to an inverse scaled sweep rate less 
than 0.15. It extends to v 2 /2tvy < 0.6 for d = 0. Increas- 
ing the population imbalance parameter d means that the 
effective density of type-2 atoms decreases, but at same 
time a type-2 atom has more chance to collide with its 
partner, the type-1 atom, because the density of type-I 
atom increases. The competition between these two ef- 
fects leads to an enhancement of heteronuclear molecular 
conversion efficiency in the nonadiabatic regime. Outside 
the nonadiabatic regime, our analytic formula overesti- 
mates the production efficiency. This deviation is mainly 
due to the difference in the range of s of Hamiltonian ([5]) 
and (|2"Tj) , respectively, i.e., it is [—1,1 — 2d] in the het- 



eronuclear case and but [—1,1] in the homonuclear case. 
This complicates the slowly sweeping case when we use 
Eq.© as the starting point of our iteration scheme. 

In the nonadiabatic regime, the conversion efficiency 
of the heteronuclear molecule can be approximated to 

. Defining the 1/e ramp rate B 1 j e 



(l+d)87r 2 H\a bg AB\n 
mB 



as that at B\/ e , T e i = 1/e, then, the ratio n/B 1 / e = 
(i+d)8eJh \ a bg AB\ is predicted to be independent of the 
density but inversely proportional to the imbalance pa- 
rameter. 

Experimentally, our theory may apply to the 85 Rb- 
87 Rb system. In Ref. |17| . the heteronuclear molecules of 
85 Rb- 87 Rb have been produced using the Feshbach reso- 
nance technique, where one BEC and a thermal gas of the 
second species are used. The main experimental parame- 
ters are a bg = 240a , AB = 4.9G, n = 1 x 10 14 cm" 3 . We 
then can calculate that the dimensionless inverse sweep- 
ing rate v 2 /(2hj) = 0.57 corresponds to the real sweep 
rate of a practical magnetic field B = 0.1G//is. How- 
ever, to apply our theory, we suggest that the experiment 
should be performed under low atom cloud densities such 
as n ~ 10 11 cto -3 with both species prepared as thermal 
gas of a few tens of nK. 



VI. CONCLUSION 

In conclusion, we have investigated the dynamics of 
the Feshbach formation of the molecules in bosonic sys- 
tems and show that the many-body effects greatly mod- 
ify the picture of two-body molecular production. With 
the help of a mean-field classical Hamiltonian, we de- 
rive an analytic expression for the conversion efficiency 
and explain the discrepancy between the prediction of 
the Landau-Zener formula and the experimental data on 
85 Kb. Our theory solves a puzzle on the formation of 
Feshbach molecules and gives some predictions on the 
formation of heteronuclear diatomic molecules such as 
85 Rb- 87 Rb. 
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